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The peculiar nature of topological surface states, such as absence of backscattering, weak anti-localization, 
and quantum anomalous Hall effect, has been demonstrated mainly in bulk and film of topological insulator 
(TI), using surface sensitive probes and bulk transport probes. However, it is equally important and 
experimentally challenging to confine massless Dirac fermions with nano-steps on TI surfaces. This 
potential structure has similar ground with linearly-dispersed photons in Fabry-Perot resonators, while 
reserving fundamental differences from well-studied Fabry-Perot resonators and quantum corrals on noble 
metal surfaces. In this paper, we study the massless Dirac fermions confined within steps along the x (r~K) 
or y(T-M) direction on the TI surface, and the Fabry-Perot-like resonances in the electronic local density of 
states (LDOS) between the steps are found. Due to the remarkable warping effect in the topological surface 
states, the LDOS confined in the step-well running along T-M direction exhibit anisotropic resonance 
patterns as compared to those in the step-well along T-K direction. The transmittance properties and spin 
orientation of Dirac fermion in both cases are also anisotropic in the presence of warping effect. 



Recently, the discovery of both two-dimensional and three-dimensional topological insulators (TIs) 1 ' 2 has 
attracted enormous attentions. Unlike the conventional two-dimensional electron states, the spin-helical 
surface states of three-dimensional TIs, which are protected by the time-reversal symmetry and consist of 
an odd number of spin-helical Dirac cones, are characterized by the gapless Dirac Hamiltonian. Owing to the 
chiral nature of the quasiparticles in TI surface, many unusual effects have been observed. Among them, the 
quantum anomalous Hall effect has been theoretically 3 predicted and experimentally 4 observed in a thin films of 
chromium- doped TI material (Bi,Sb) 2 Te 3 , which may open a door of the application of TI materials in the field of 
the low-power-consumption electronics. A robust response of weak antilocalization 5 " 7 has been tested in the 
transport studies, and the suppression of the backscattering has also been confirmed by examining the scattering 
from impurities or step edges 8 " 18 in previous scanning tunneling microscopy (STM) reports. Particularly, because 
of the warping effect 19 in topological surface states, the step scattering induced Friedel oscillations in the local 
density of states (LDOS) exhibit quite unique decaying power law 1112 as compared to that in noble metals 20 " 23 and 
graphene 24 . Furthermore, with the advanced progress of STM technologies, it is possible to study the optics- 
analogous properties of the quantum interference as well as quantum confinement of the conventional two- 
dimensional electron gas on noble metals 25 ' 26 . Recently, many efforts have also been devoted to investigate the 
Dirac electron analogue of optical phenomena, such as negative refraction, Goos-Hanchen effect, beam collima- 
tion, and Fabry-Perot resonances/interferences, in graphene systems 27 " 34 and TI materials 10 ' 35 " 37 . For instance, 
Fabry-Perot type conductance oscillation patterns of Dirac fermions in monolayer graphene p-n junctions 30,34 
and giant conductance oscillations in ballistic trilayer graphene Fabry-Perot interferometers 32 have been suc- 
cessfully observed in recent quantum transport experiments. More recently, in contrast to ordinary surface states 
in common metals, recent STM experiment performed on TI material Sb(lll) surface suggests that the Dirac 
fermions transmit the step barriers with a high probability, and the Fabry-Perot resonance of the topological 
surface states has been observed 10 . Although it is highly challenging in experiments to confine the massless Dirac 
fermions in the nanoscale quantum corrals 38 constructed by adatoms or Fabry-Perot resonators formed by 
straight parallel step defects on TI surface, the latest success has been achieved on Bi 2 Te 3 (lll)-based samples 39 . 

Because of its importance both from basic point of interest and to TI -based quantum device applications, in the 
present paper we address this issue by presenting an attempt at the theoretical evaluation of the Fabry-Perot 
resonance problem of the massless Dirac electrons on the TI surface in the presence of double symmetric steps. 
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We show that because of the strong warping effect in the topological 
surface states, the electronic LDOS confined in the step-well running 
along y (r~M) direction and those in the step-well running along 
x (r-K) direction exhibit anisotropic Fabry-Perot resonance images. 
The resonant transmission properties as well as the spin orientation 
of Dirac fermion in both cases are also influenced remarkably by the 
warping effect. These findings could be confirmed by STM 
measurements. 

Results 

We start from the effective Dirac Hamiltonian for the TI surface 
which is expressed as 1119 ' 40 

H 0 (*) = v/ {a x ky - o y k x ) + ^ (k\ +k 3 _)a z , (1) 

where (—240 meV'nm) is the Fermi velocity, X (—250 meV'nm 3 
for Bi 2 Te 3 ) is the warping parameter 9 ' 1119 , & = {o x , o y , o z } are Pauli 
matrices acting on spin space, and k± = k x ± ik y . For simplicity, here 
we ignore the spin-independent k 2 /2m term in the effective 
Hamiltonian since it just results in particle-hole asymmetry but 
affects the shape of Fermi surface little. Therefore, the anisotropic 
Fabry-Perot resonant states discussed following will not be affected 
when the particle-hole asymmetry is absent. The eigenfunction of H 0 

is given by \j/ s (k, r) = (<j) u , si(j)2, s e ie ) r e ikr , where e 10 = x y , and 

K 

s = ±1 corresponds to the upper and lower band dispersion 
e s (k) = syj (vfk) 2 + k 2 k 6 cos 2 (30), which is schematically shown in 

k! 



d — sds 
Id 



where 



k' x , (k'\ + k f3 _ ) jl\ with k' = kJI/vf = kl 



Fig. 1(c). Here, (j) ls = 

d = {^1,^2^3} — {k'y, - 

Usually, the straight and parallel steps could be naturally formed 
in the progress of molecular-beam epitaxy growth, which will result 
in scattering of quasiparticles on the surface of the sample. Hereafter 
we consider the problem of massless Dirac electrons confined in a 
pair of symmetric steps apart la along the x or y direction on the 
surface of a three-dimensional TI, which is schematically shown in 
Fig. 1(a). The step-edge potential is also illustrated in Fig. 1(b), which 
is assumed as one-dimensional (5-type. 

Steps along the y direction. In what follows, let us consider the 
double nanoscale parallel steps along the y direction on TI surface 
with the scattering potential described as (5 -type potential 

barriers 15 ' 16 ' 20 " 22 



U(x) = Vh [S (x + a) + S (x - a)] , 



(2) 



where h is the hight of step and the step edges are located at x 0 = ±a. 
At this stage, we would like to point out that it is easy, on one side, to 



(a) 



(b) 



U 



m 




Figure 1 | Two symmetric steps on a TI surface with strong warping 
effect, (a) Schematics of the symmetric steps and (b) double (5-type 
barriers, (c) Constant energy contour of e 5 (k) for Bi 2 Te 3 TI material with 
the fermi velocity Vf = 240 meV'nm, and the warping parameter X = 
250 meV'nm 3 . 



treat the boundary conditions for a pair of (5 -type potential barriers, 
while for other types of barriers (such as rectangular, semielliptic, 
and Gaussian profile barriers) the eigenequations of the boundary 
condition become too tedious and complicated to obtain the 
analytical results. On the other side, the key physical properties in 
the step scattering will not be lost by using (5-type barriers herein. 
Therefore, in this work we choose to employ barrier model (2) for 
simplicity and clarification. 

In the scattering process, we suppose an incident electron plane 
wave from one side [region I in Fig. 1(a)] with momentum k and 
energy £ s will be reflected back to the same side or transmits into the 
step- well (region II), and then the transmitted part will be reflected 
between the both steps or transmits further into the other side 
(region III). Thus, the wave functions could be written as 



^(r) = *A 5 (*,r) + n^ s (k f ,r) , (x<-a) 
^n(r)=A^ s (Kr) + B^ s (k f ,r) (\x\<a) 
il/ in (r)=Txl/ s (k,r), (x>a) 



(3) 



Here, we have assumed that k y is a good quantum number and the 
energy is conserved in the electron scattering process, naively, the 
reflected and transmitted waves will be characterized by kf = (~k x , 
k y ) = {k y Of) and k = (k x , k y ) = (k, 0) with 6 f = n - 6. Such an 
assumption is an analogy between electron transport and light pro- 
pagation since it could be understood in terms of phenomena like 
reflection and transmission 31 . This results in that 0i, 5 (— k x , k y ) = 
<t>i- s (k x , k y ) and (j)2,s(~k x , k y ) = (/>2- s (kx> ky)- H an d T denote the 
reflection amplitude and tunneling amplitude, respectively. By con- 
sidering the boundary conditions of wave functions at the step edges 
Xq = ±a, we can obtain the amplitudes A, B, 71, and T, seeing the 
details of the derivation in the section of METHODS. 

It is also easy by using Eq. (3) to get the analytical expressions of 
spin orientation (<r) = (i// (r)\o\\jj (r)) of Dirac fermion on TI surface 
in the presence of Fabry-Perot resonator, which are shown in the part 
of METHODS. Because of the scattering and interference induced by 
the double steps, the spin texture of the bands becomes spatially 
dependent, and is modified via the amplitudes A, B, 7Z, and T, which 
could be clearly found from Eqs. (12-14). As an example, in Fig. 2(a) 
we exhibit the spatial distribution of spin orientation for a scattered 
Dirac fermion with incident angle 6 = n/3. It is obvious from 
Fig. 2(a) that (v x ,y, z ) x <-a and (v x ,y, z )\ x \< a in incoming and step-well 
regions [i.e., region I and II shown in Fig. 2(a)] periodically oscillate 
with space coordinate x while ((T x ,y,z)x>a m out-coming region 
[region III in Fig. 2(c)] are independent on x, which are consistent 
with the analytical results in Eqs. (12-14). Furthermore, one can also 
observe from Fig. 2(a) that the spin orientations of ((T Xy y yZ ) are dis- 
continuous at the steps x 0 = ±a indicating sudden change of the 
direction of spin, which is similar to the case of a single step (not 
shown here for briefness). The numerical results of spin orientation 
corresponding to Eqs. (12-14) as a function of the direction of in- 
plane wave-vector are shown in Figs. 2(c-e), and for comparison, the 

spin orientation for the free Dirac fermion on a clean surface of TI 
5 

(a) = - (k' sin0, — k! cos0,d 3 ) is also shown in Fig. 2(b). In Figs. 2(b- 
d 

e), the black arrows denote the xy-plane component (<T||) = {o x ) + 
i{<j y ), the color background denotes the z component (a z ), and the 
wine curves denote the contours of constant energy. We can see from 
Fig. 2(b) that for the free Dirac fermion, the xy-plane component 
(<7||) rotates clockwise (s = + 1), and the threefold symmetry of (a z ) 
becomes obvious with increasing the Fermi energy due to the warp- 
ing effect. Comparing with the case of a free Dirac fermion 
[Fig. 2(b)], due to the quantum interference, the spin orientation 
of the scattered Dirac fermion located in regions I and II is affected 
so strongly that on one hand, not only the amplitude of (<7||) is 
changed remarkably but also the clockwise rotation of (<T||) is blured, 
on the other hand, the threefold symmetry of (a z ) is broken, see 
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Figure 2 | Spin orientations of Dirac fermions on the surface of TI. (a) The spatial distribution of spin orientation for Dirac fermions scattered 
from the double steps, where the incident angle 6 = n/3. (b) Spin orientation of free Dirac fermions. (c-e) Spin orientation of scattered Dirac fermions at 
Region I x = — 1.5a (c), region II x = 0.5a (d), and region III (e) on a TI surface with double steps located at ± a. In (b-e) the arrows (color backgrounds) 
denote the xy-plane component (<ti i) (^component (cr z )), and the wine curves denote the constant energy contours. Parameters are chosen as a = 20 nm, 

Vh 

s = +1, Vf= 240 meV'nm, warping effect X = 250 meV'nm 3 , and w 0 = — - = 1.5. 



Figs. 2(c) and 2(d). Whereas, for the Dirac fermion transmitted into 
region III, there are two key points that need to be noticed: i) Both the 
clockwise rotation of (<7||) and the threefold symmetry of (a z ) are 
exact; ii) The amplitude of {o x ^ z ) x>a is proportional to the transmit - 
tance |T| 2 [Eq. (14)], and thus (<j z ) exhibits sharp peaks in Fig. 2(e) 
which is corresponding to resonant tunneling (see the following 
discussion). Moreover, in the absence of warping effect of surface 
bands (i.e., d 3 = 0 and k' = d), the spin orientation (a) could be 
further simplified (see the formulae below Eq. (15) in the METHODS 
section). Especially, we find in this case that i) (<7 y ) will not oscillate in 
real-space since it becomes independent on x, and ii) the spin of the 
Dirac fermion in region III is still locked in the surface plane since 

(°z)x>a = 0. 

Now let us turn to discuss the transmission property of Dirac 
fermions in the Fabry-Perot resonator on TI surface, which is also 
irradiative for the explanation of the spin orientation discussed 
above. Differing from the single step case 15 " 17 , there will appear res- 
onant tunneling in two-step or multi-step scattering. We could ana- 
lytically determine from Eq. (15), that the resonance transmission 

mi 

satisfies the condition oik n v = — with n an integer number. In other 

words, for the incoming state with wave vector k = (k, 0), the res- 
onance transmission occurs when the incident angle satisfies the 

relationship that cos 9 = — - , and the corresponding 



( nn 2 



resonance level is expressed as e n;5 (k",k y ) = svj y ^— + kjj in the 

absence of warping effect in the topological surface state. It is easy to 
check that the reflection coefficient and tunneling amplitudes satisfy 



the relation of |7£| 2 + |T| 2 =.R+T=1. Furthermore, we should 
point out that if the plane wave incident from the right side of steps, 
i.e., from region III, one could also obtain the similar expressions of 
W and T, and \K'\ 2 + \T\ 2 = \U\ 2 + |T'| 2 = 1 is satisfied, which is 
the same as that of a single step case. Similar to the case of the optical 
Fabry-Perot resonator, 1Z and T are also determined by the distance 
between the two mirrors and the incident angle of light, and if the 
absorption is neglected |7£| 2 + |T| 2 = 1 is required because of the 
conservation of energy. However, differing from the optical case, 
1Z and T of Dirac fermions in the present system are strongly 
dependent on the nature of bands (warping effect) and the height 
of barriers seeing the following discussions and formulae in the 
METHODS section. 

The typical results of reflectance and transmission strengths as 
functions of incident angle 0 (the angle of incident wave vector k 
to the normal direction of steps) are shown in Fig. 3. One can find the 
resonant tunneling T = 1 (correspondingly, the backscattering is 
forbidden) at normal incidence (6 = 0) and other special values of 
6 that satisfy the resonance conditions, see the red curves in Fig. 3. 
Another interesting phenomenon should be noticed that although 
the spin orientation of Dirac fermion rotates in the step -well region 
due to the quantum interference, when the resonant tunneling 
occurs, its orientation in the out-coming region will turn back to 
that in the incoming region, i.e., {a) x>a = {a) x<a = 

5 

-(/c'sin0, — k' cos6,d 3 ) when T = 1 and R = 0, see Eqs. (12) and 
d 

(14). In other words, the spin of transmitted Dirac fermion points 
along the same direction as that for the incident Dirac fermion when 
it resonantly tunnels through single- or multiple-step. It can also be 
found from Fig. 3 that with increasing the energy of incident electron, 
the resonant tunneling becomes more remarkable. By a comparison 
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Figure 3 | Transmittance (red) and reflectance (blue) strengths of Dirac fermions on the surface of TI with two symmetric steps grown along 

the /direction, (a-b) Transmittance Tand reflectance R vs. the incident angle 6 without warping effect, (c) Tand R with warping effect. The momentum 

is chosen as k = 0.2 nm" 1 in (a), and k = 1.3 nm" 1 in (b) and (c). Other parameters are the same as those in Fig. 2. 



of Fig. 3(b) and Fig. 3(c), one can see that the curves of reflectance 
and transmittance strengths are reformed due to the influence of 
strong warping effect, but the incident angles corresponding to the 
resonant tunneling are unaltered. 

In order to determine the discrete resonance energy levels con- 
fined in the step-well, we plot the transmittance strength versus the 
wave vector components k x and k y without [Fig. 4(a)] and with 
[Fig. 4(c)] warping effect. Correspondingly, the colored background 
in Figs. 4(b) and 4(d) are the transmittance strength versus the 
energy s of electron and the wave vector component k y . With increas- 
ing energy, it is clear that the warping effect becomes strong and 



obvious. We can observe the resonant tunneling peaks correspond- 
ing to the quantized k x , and from these numerical results of T the 

YlTl 

resonance condition can be also concluded as k i l= — , which is 

x 2a 

consistent with the analytical results obtained from the expression 
for T in the absence of warping term [see Eq. (15)]. Substituting this 
k n x back into the band dispersion s s (k), one can obtain the discrete 
resonance energy levels s n , s (k" 9 ky) confined in the step-well, which 
are also shown in Figs. 4(b) and 4(d) by the black dotted curves. It is 
clear that the resonance e n?5 (k x ,k y ) are perfectly consistent with the 
resonant tunneling peaks. 




-1.0 -0.5 0.0 0.5 
k(nm- 1 ) 



k(nm ) 



Figure 4 | Resonance transmission of Dirac fermions on the surface of TI with two symmetric steps grown along the y direction, (a-b) Transmittance 
contour plot versus k x and k y , and versus k y and e for the case without warping effect, (c-d) Transmittance for the case with warping effect. 
The black dotted curves in (b) and (d) are for the resonance energy 8 n>s . The parameters are the same as those in Fig. 2. 
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After obtaining the resonance energy levels, one can easily get the 
LDOS between the double steps, which is expressed as 



k»s 



dk 



(4) 



In order to observe the Fabry-Perot resonant states more clearly, we 
subtract the contributions of the averaged background p° u and only 
consider the change in the LDOS Ap n (r, e) due to the step scattering 
in following calculations. 

Figure 5 shows the energy and position dependence of the LDOS 
confined in a symmetric step-well along they (r~M) direction on TI 
surface. One can clearly see the obvious quantum confinement, and 
the LDOS images are symmetric because our Fabry-Perot step-well 
model is symmetric in the calculations herein. The number of res- 
onant peaks increases with increasing energy, which is consistent 
with the resonance properties of transmittance shown above. In 
our calculations, the ground state confined in the Fabry-Perot res- 
onator is located at 40 meV, while the first excited state with one 
node lies at 55 meV. When comparing Fig. 5(a) with Fig. 5(b), it is 
surprising to find that there exists a strong enhancement in the LDOS 
in the energy region from 120 meV to 200 meV when the warping 
effect is taken into account in the calculations (X 7^ 0), while at other 
energy regions the resonance modes are similar to the case in the 
absence of warping effect (X = 0). The differences in peak intensity 
and resonance period could be more clearly observed from the curves 
in Figs. 5(c)-(e). In addition, we have also checked other step-wells 
with different widths, the enhancement of the LDOS around the 
energy region of 120 ~ 200 meV is still remarkable (not shown here). 
This enhancement is closely related to the occurrence of the flat-band 
crossover when the Aylinear positive term, which is the dominate 
term at low energy, is superposed onto the Ay square negative term 
which is dominate at high energy [see Fig. 3(d)]. This bunches the 
electronic states and therefore results in the enhancement of the 
LDOS in the observed energy region. One could also qualitatively 
understand the enhancement of LDOS in the energy region of 120 ~ 
200 meV from the view point of the approximation of stationary 
phase points (SPPs) 11 . The LDOS is dominated by the SPPs [(k x , 
k y ) and (—k x , k y )] along the direction of T-K on the constant energy 
contour since most of the contributions to the integral in Eq. (4) are 



from the vicinity of these points. Therefore, the integral can be bro- 
ken up into regions between stationary phase points that 



Ap n (x,e)~Re 



AB*- 



M 1- 



0 2id\ 



2dn 



dk v 



where 



/. 



means that the integration of k y is around SPPs. Taking into account 
the warping effect, the separation between the resonance levels at 
energy region of 120 ~ 200 meV (n = 5, 6, 7, 8) is compressed, see 
the black dotted curves in Fig. 4(d). Thus, the LDOS will be 
enhanced. While with further increasing of energy, the separation 
between resonance levels becomes larger and larger, which results in 
the decrease of LDOS [see Fig. 5(b) and the black curves in Figs. 5(c) 
and 5(d)] . Additional q- vectors (q = kf— k t ) apart from the SPPs will 

blur the LDOS because of the rapidly varying phase factor 

in the integral of Eq. (4), which has little contribution to the total 

LDOS. 

The numerical results of LDOS outside the step-well Ap 7 (r, s) are 
plotted in Fig. 6. We can see from Fig. 6(a) that resonant tunneling of 
the Dirac fermions is obvious in the contour plot of energy- resolved 
modulations of the LDOS, which is consistent with the above dis- 
cussions of the transmittance. Different from the single step case, the 
quantum interferences due to the multi-step scattering is clear in the 
present Fabry-Perot resonator, and the warping effect is also remark- 
able [comparing the black curves with the blue curves in Figs. 6(b) 
and 6(c)] with the increase of energy. In addition, the character of the 
Friedel-type decaying oscillations in the present case is similar to the 
single step case. Our theoretical simulations are qualitatively consist- 
ent with the STM experimental observations of the Fabry-Perot 
resonance of the topological surface states on the Sb(lll) surface 10 , 
and the theoretical methods used in this work should be useful for the 
explanations of the related STM experiments. One should also notice 
that because k y keeps unchanged in the scattering process, the LDOS 
oscillates only in the perpendicular direction of steps, which is the 
same as the spin orientation. 

Steps along the x direction. Now we switch gears and consider the 
case of steps along the x (T-K) direction. Due to the hexagonal 
warping term in the Hamiltonian (1), the x (T-K) and y (r-M) 
directions are inequivalent, which could be seen from the constant 
energy contour shown in Fig. 1(c). Moreover, the warping term in 



tu 400 u i 

(b) 




;(c) t 


x=0.0 









100 200 300 400 
c(meV) 
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■ -0.04__ 

-20 -10 0 10 20 -20 -10 0 10 
x(nm) x(nm) 



Figure 5 | The change in the LDOS of Dirac electrons confined in the double symmetric steps apart along the y direction on TI surface, (a) The 

change in the LDOS Ap u (x, e) confined in the double symmetric d-type steps apart 2a = 40 nm on TI surface without warping effect, (b) Ap n (x, co) 
confined in the steps on TI surface with warping effect, (c) Bias dependent LDOS spectra at the center of step-well x = 0. (d-e) Spatially dependent LDOS 
spectra within finite bias e = 150 meV and e = 350 meV, respectively. The blue (black) curves are for the case of without (with) warping effect. All the 
other parameters are the same as those in Fig. 2. 
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Figure 6 | The change in the LDOS outside of the double steps along the y direction on TI surface with warping effect, (a) The change in the 
LDOS Api (x, e) versus x and e. (b) and (c) spatially dependent Apj (x, e) spectra within finite bias e = 150 meV and e = 350 meV, respectively. The blue 
(black) curves are for the case of without (with) warping effect. 



Hamiltonian (1) is represented as i^{d x — d x d 2 ^j o z in the real space, 

therefore, it is naturally expected that different orientation of steps 
leads to different Fabry-Perot resonant features in the LDOS and in 
the reflection/transmittance properties. For the steps along the x {T- 
K) direction, the wave function in each region should be written as 



^{r) = ^ s {k,r) + n^ s (k f ,r), 



(y<-a 
(\y\<a) 
{y>a) 



(5) 



In this case, the reflected wave vectors will be characterized by kf = 
(k x , -k y ) = (k, -0), which results in (j) hs (k x , -k y ) = (fri, s (k x , k y ) and 
fyijkx* ~k y ) = fyijjkx* k y ). Following the same method used above to 
consider the boundary conditions aty 0 = —ci, one could also easily 
get the coefficients A, B, 71, and T [see Eq. (16) shown in the section 
of METHODS]. 

The reflectance and transmittance strengths for this case are 
shown in Fig. 7. Similar to the case of steps along)/ (r~M) direction, 
we can observe the resonant tunneling herein, and the resonance 

YITC 

condition could be concluded as k n v = — . However, the transmit - 

y 2a 

tance mirages and the quantized resonance energy £ n ,s(^ x ,ky^j are 

different from those for the steps along y (r~M) direction, which 
could be clearly observed by comparing Fig. 7 with Fig. 4. Here, the 
quantized energy dispersion is a result of the superposition between 
the /^-linear positive term and the /c x -cubic positive term. Therefore, 
the aforementioned flatband crossover cannot occur in the present 
case, as verified in Fig. 7. Notice that 0 in Fig. 7 is defined in the same 
way as in Fig. 4, i.e., the angle of incident wave vector to the positive x 
axis. Besides, one could derive the spin orientation in the present 
case, and the anisotropy could also be found in the spin orientation 
since it is related to the coefficients A, B, 7Z, and T (not shown for 
briefness). 

The LDOS results for this case are illustrated in Fig. 8. When the 
warping effect is taken into account, the Fabry- Perot-like resonance 
of LDOS confined in the step-well along the x (r~K) direction is 
quite different from the case of steps along the y (r~M) direction, 
which indicates the anisotropic step -scattering character on the TI 
surface. In this case, the LDOS is dominated by the SPPs along T-M 
direction on the constant energy contour. Especially, the Fabry-Perot 
resonance is enhanced [see Fig. 8(b)] around 300 ~ 350 meV in this 
case since separation between the resonance levels at this energy 
region is compressed [see the black dotted curves in Fig. 7(c)] due 



to the warping effect. These interesting anisotropic Fabry-Perot res- 
onance states confined within nano- steps on the TI surface can not 
be observed on conventional metal surface or graphene. 

Discussion 

Previous STM experiment 9 and theory 11 about the scattering of sur- 
face states off single step edges in TI materials Bi 2 Se 3 and Bi 2 Te 3 in 
which the warping effect is strong suggested that the power law decay 
of the standing waves could be modulated with increasing the Fermi 
level due to the strong warping effect. Based on the above results, we 
would like to clarify that the experimental relevance of the aniso- 
tropic Fabry-Perot-like resonances states on TI surface will also 
strongly depend on whether the Fermi level lies close the Dirac point 
or lies far from the Dirac point so that transport at Fermi level can get 
strongly effect by warping. The modulation of Fermi level could be 
easily realized at present experimental conditions by appropriate 
doping. Chen et al. 41 experimentally suggested that the Fermi level 




-0.5 0.0 , 0.5 
k^nm ) 

Figure 7 | Transmittance and reflectance strengths of Dirac fermions on 
the surface of TI with two symmetric steps grown along the x direction. 

(a) Transmittance T and reflectance R versus the incident angle 6. 

(b) Transmittance contourplot vs k x and ky, and (c) transmittance 
contourplot vs k x and e (c). The warping effect is taken into account in 
calculations, k = 1.3 nm" 1 in (a). The black dotted curves in (c) are for the 
resonance energy e n fc . 
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of Bi 2 Te 3 can be tuned to intersect only the surface states with strong 
warping effect, and a full energy gap for the bulk states was clearly 
observed in the angle-resolved photoemission spectroscopy experi- 
ment. This indicates that as the Fermi level goes away from the Dirac 
point, although the angle between spin and momentum is no longer 
locked in the surface plane since an additional component of out of 
plane orientation of spin is raised by the warping effect, the bulk 
disorder may play a weak role in the scattering process of topological 
surface states. Thereby, we have not considered the role of the bulk 
disorder in our calculations. 

Furthermore, on one side, similar to the case of scattering of sur- 
face states caused by single-step defect on the surface of TI 15 " 17 , 
anisotropic transmittance and decay of LDOS outside of the nano- 
step well has also been observed in our studies due to the strong 
warping effect; On the other side, however, apart from various details 
of the model, the major distinctive features of our study are (i) the 
resonance and interference behaviors of Dirac fermions transmitting 
double step barriers as well as the corresponding resonance energy 
levels [Figs. 3, 4, 7], and (ii) the anisotropic Fabry-Perot resonant 
states confined within nano-step wells grown along different high 
symmetry directions on TI surface in the presence of strong warping 
effect in surface states. Recently, an STM measurement showed that 
the Fabry-Perot resonance of the topological surface states on 
Sb(lll) surface with multiple asymmetric steps has been successfully 
observed 10 . The high-quality Fabry-Perot interferences of Dirac fer- 
mions in graphene-based devices have also been achieved in recent 
quantum transport experiments 30 ' 32 ' 34 . More recently, both the 
quantum confinement of massless Dirac fermions in nanoscale 
quantum corrals surrounded by Bi-bilayer on Bi 2 Te 3 surface 39 and 
the anisotropic scattering of surface state electrons at a point defect 
on Bi(lll) surface 42 have been experimentally detected. These 
experimental observations indicate that our methods should be help- 
ful for understanding these experimental findings, and we expect that 
the anisotropic Fabry-Perot resonance states found in this work 
could be confirmed in future STM experiments. 

In summary, the Fabry-Perot-like resonances of massless Dirac 
electron confined within symmetric steps along the x (r-K) or y (T- 
M) directions on TI surface have been studied. The analytical expres- 
sions for the spin orientations, reflectance, and transmission were 
obtained in the presence of the warping effect. We found that the spin 
orientation of the Dirac fermion rotates in the step-well region due to 
the quantum interference, but the transmitted Dirac fermion points 
along the same direction of the incident fermionic spin when it 



resonantly tunnels through the step -well. Particularly, because of 
the strong warping effect in the topological surface state, the spin 
is no longer locked in the surface plane, and the electronic LDOS 
image confined in the step-well running along T-K direction is 
remarkably different from that in the step -well running along T- 
M direction. The resonant transmittance properties as well as the 
spin orientation of topological surface states in both cases are also 
different due to the warping effect. The formula obtained here can be 
extended to symmetric/asymmetric multi-step cases as well as spin 
polarized cases. Also, our anisotropic results, which can not be 
observed in the conventional metal surface and graphene system, 
may provide new insights in constructing a TI -based Dirac electron 
Fabry-Perot wavevector filter, which is analogous to the optical 
Fabry-Perot filter. While differing from the optical Fabry-Perot cav- 
ity, which could be used to store photons in a gravitational wave 
detection, the Fabry-Perot cavity discussed here may be useful for 
confining and storing Dirac electrons. 

Methods 

We derive an equation of continuation of wave functions at boundaries to get the 
coefficients A, B, 7Z, and T for the case of the steps grown along the y (IVM) 
direction, which is used to calculate the LDOS in our system. Within the Eq. (3), one 
could get the total wavefunction as 

ij/ tot = il/ I @(-a-x) + il/ n [@(a + x)-®(x-a)} + il/ III ®(x-a), (6) 

where 0(z) is the step function. Now we consider the property of wavefunction at x 0 
= ± a respectively. Substituting Eq. (6) into the Hamiltonian ( 1 ) in real space and plus 
the 3 -type potential, we have 

[U (d\ - 3d x d^j o z - iv f d y o x + iv f d x a y ^ \J/ tot (x,y) 

= {8-Vd{S(x + a) + d(x- a)}}ijy t0t (x,y) . 

Integrating both sides of the Eq. (7) near the points x 0 = ±a, 

lim [ dx \ik (d\ - 3d x d 2 y ) o z - iv f d y o x + iv f d x o y ] \jj tot (x,y) 

8 

r±a+e V ' 

= lim dx[£ — Vdd(x + a)]\j/ tot (x,y), 

e ^°J±a-e 

we have the following algebraic equations 

' aQ+a l K + a 2 A + a 3 B = 0 
b Q + b l n + b 2 A + hB = 0 

ciT + c 2 A + c 3 B = 0 ' 1 1 

k /iT+m+/ 3 £=o 

where the parameters are expressed as 
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a 0 = 



b Q = [{u-w)s^ s e ie -<t> u ]e 



a x = [( w + u )(l) h _ s + s(l) 2 _ s e- i6 ]e ikxa , b x = [{w-u)s<j) 2 _ s e- i6 -<j) h _^ kxa , 
a 2 =[(u-w)(l) u + s(j> 2 / e ]e- ik * a , b 2 = [{w + u)s(l) 2>s e ie + <l) hs ]e- ikxa , 

a 3 =[(u-w)4> h - s -s^_ s e- i6 ]e ik * a , b 3 = [- (w+ u)s^_ s e~ lG + ^_ s ] e^\ 

Cl = [{u-w^+s^y^, u = [( w+M )5^+^j^, (io) 

c 2 = [{w+u)^ u -s^/y k ^ a , f 2 =[(u-w)s^/ e -^2e^ a , 

c 3 = [(w+u)(\)^_ s + s(l) 2 _ s e- ie ]e- ikx \ f 3 = [(w- u)s(j) 2 _ s e- i0 -^i_J e~ ikxCl . 

From Eq. (9), one can get the coefficients A, B, 7Z, and T. In the presence of the 
warping effect, the coefficients are easily written as 



M M M M 



(11) 



where M = (a 2 b x - a x b 2 ) (c 3 f x - c x f 3 ) - (a 3 b x - a x b 3 ) (c 2 f x - c x f 2 ),y = (a 0 b 3 - a 3 b 0 ) 
(cxfi ~ c/i) + (a 2 b 0 - a 0 b 2 ) (cj 3 - c 3 f x ), cc = (a 0 b x - a x b 0 ) (c/ 3 - c 3 f x ), /3 = (a x b 0 - 

Vh 

«o&i) {C]f 2 - c/i), and t = (a x b 0 - a 0 b x ) (c^f 3 - c 3 f 2 ). Here, u = —— = -u 0 i and 

ZlVf 

w= (k 2 x — 3k 2 ^J I 2 are dimensionless variables. We should point out that the prop- 

pa + e 1 

erties of Dirac <5-function that lim e ^ 0 dxS(x — a)f(x) = - \f (a + ) +f(a )], and 

Ja—e 2 

have been used in the derivations of Eq. (9). 

After obtaining the wavefunctions, we could easily get the analytical expressions of 
spin orientation {a) = (\J/(r)\(7\\J/(r)) of Dirac fermions on TI surface in the presence 
of Fabry-Perot resonator. In the region I (incoming region) in Fig. 1(a), we have 



J dxf(x)S^ (x) = - [ dx 3 -^- d^' l \x) 



sk' 



1 + \n\ 2 )sm8 + 2slm[n.*e i{2k ' x+9 >} , 
^|T| 2 cos0+2^Re[RV< 2 ^ +, '>], 



(12) 



(*z) x< -» = S ±\T\ 2 + k jHK'e 2ik <*(e™ + l)}. 
Similarly, in the region II (step-well region), we have 

{^> w<a =^(|^| 2 + |S| 2 )sin0 + 25lm[^e^+")], 
W) l x|<«, = y(|S| 2 -M| 2 )cos9 + 2|Re^«^")] 
(^) lx]<a = ^Re[AB*e 2 >"-*(e 2 * + l)} 

+ ^ { \A\ 2 - \B\ 2 -Re[AB'e 2ik '*{e 2i0 - 1)] }, 
and in the region III (out-coming region), we have 



(13) 



<»r>,>. 



|T| 2 cos0, 



(14) 



In particular, if the warping effect is ignored (k = 0) in the derivations, the coefficients 
should be reduced to the following simpler form that 

s(l-w 2 )cos0[2w + s(l + w 2 )cos 9] 
Mo ' 



se i ^ a - 6 )u{l-u 2 ){e m -l) 



11 = 



2M 0 

-l)[ r _g2^a + r+e - 



(15) 



M 0 



B + l)\l-u 2 



4M 0 

where M 0 = [2w + s(l+w 2 )cos0] 2 -2e 4iM w 2 [l-cos(20)],andr ± =5(1 + u 2 ) cos 9 
± 2. Correspondingly, the spin orientation without considering the warping effect 
could be simplified, and we have in region I 

(( 7jc ) x< _ fl =s(l + |^| 2 )sin0 + 2sIm[^*e /(2 ^ + 0) ], (a y ) x< _ a = -s\T\ 2 cos9, and 
(o z ) x< - a = Re [K*e 2ikxX (e 2i0 + 1)] , in region II 
K> W <a = Kl^| 2 + l^| 2 )sin0 + 25lm[^*e^ +0 )], 



(a,) w<fl =5(|^| 2 -|^| 2 )cos0,and (a z ) lx{<a = Re[AB* e> ik * x (e 2ie + 1)], and in 
region III (o x ) x>a =s\T\ 2 sm9 y {cy) x>a = -s|T| 2 cos0, and ((J z ) x>a = 0. 

Furthermore, for the case of the steps grown along the x (r-K) direction, with the 
same method above to treat the wave functions (5), we could also get the coefficients 

y - oc - B - t 

1Z = — , A = — , B= — , and T= — , where M, a, B, y, and x have the same form 

M M M M 
shown above but with the parameters expressed as 



a 0 = 


[ - (w+u)(\) X s - s4^/ 0 ] e~ ik y a , 


b 0 = 


[~<i>i,s + ( w ~ u ) si hs eW ] e~ ikya , 


a x = 


[{w-u)^-si^ s e- i6 Y k y\ 


h = 


[-(w+ W )5^ 2/ - !0 -^J e V, 


a 2 = 


[(w-u)^ s + si(j> 2 ^]e- ik y\ 


b 2 = 


[-(w+u)si<j> 2 / B + <j> hs ]e- ik > a , 


a 3 = 


[-(w+u)(j) x ^ + si^ s e- ie y k y a , 


b 3 = 


[(w-u)si^ s e-' d + ^ s ]e lk y a , (16) 


C\ = 


[{w-2u)^ s + si^ 2 / B y k y\ 


h = \ 


- (w+2u)si<j) 2 je- ie + <^ 1)S ] e ik y\ 


c 2 = 


[-(w+u)$ hs -si^ s e l6 ]e ik y a , 


h=[ 


\w-2u)si(j) 2 y e -(j) x ^ e ik y\ 


C3 = 


[(w-2u)(]) Xs -si(]) 2s e- id ]e- ik y a , 




- ( w + 2 u) si0 2>5 e ~ ie - (j) x J e ~ ik y a . 



After obtaining these coefficients, one could get the LDOS confined in the nanosteps 
by employing Eq. (4) in the text. 
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